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There is much latitude between the requirements of 
Schnirelman's theorem regarding the ergodicity of individual 
high-energy eigenstates of classically chaotic systems on the 
one hand, and the extreme requirements of random matrix 
theory on the other. It seems likely that some eigenstate 
statistics and long-time transport behavior bear nonrandom 
imprints of the underlying classical dynamics while simultane- 
ously obeying Schnirelman's theorem. Indeed this was shown 
earlier in the case of systems which approach classical ergod- 
icity slowly, and is also realized in the scarring of eigenstates, 
even in the h ^ limit, along unstable periodic orbits and 
their manifolds. Here we demonstrate the nonrandom charac- 
ter of eigenstates of Sinai-like systems. We show that mixing 
between channels in Sinai systems is dramatically deficient 
compared to random matrix theory predictions. The deficit 
increases as |logft| for h —^ 0, and is due to the vicinity of 
the measure zero set of orbits which never collide with the 
Sinai obstruction. Coarse graining to macroscopic scales re- 
covers the Schnirelman result. Three systems are investigated 
here: a Sinai-type billiard, a quantum map which possesses 
the essential properties of the Sinai billiard, and a unitary 
map corresponding to a quasirandom Hamiltonian. Various 
wavefunction and long-time transport statistics are defined, 
theoretically investigated, and compared to numerical data. 



I. INTRODUCTION 

In recent years, much attention has been paid to 
the structure of quantum eigenstates in systems with 
a chaotic or ergodic classical analogue. For integrable 
systems, EBK quantization provides an intuitive under- 
standing of classical-quantum correspondence, associat- 
ing quantum wavefunctions with the invariant tori of the 
underlying classical dynamics. In a classically ergodic 
system, the typical trajectory fills an entire energy hyper- 
surface at long times, and it is natural to conjecture that 
a typical high-energy eigenstate of such a system simi- 
larly has intensity distributed evenly over an entire en- 
ergy shell. Thus, Berry suggested in 1983 that an eigen- 
state of a classically ergodic system should look locally 
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like a random superposition of plane waves of fixed en- 
ergy, with momenta pointing in all possible directions [El . 
Similarly, Bohigas, Giannoni, and Schmit [^ proposed 
that the quantum properties of a classically chaotic sys- 
tem should correspond to those of random matrix theory 
(RMT). This implies that wavefunction intensity should 
be distributed over an entire energy surface, with the 
wavefunction amplitudes at distant points behaving as 
independent Gaussian variables. 

The conjecture that chaotic eigenstates obey RMT 
statistics is a statement about quantum structure at the 
scale of a single wavelength in position space (or on 
the scale of a single channel in momentum space, or 
more generally, on a mesh of cell size 0(h) in phase 
space). Rigorous results on quantum ergodicity, how- 
ever, mostly address structure on classically large scales, 
in the limit where h becomes small compared to the 
phase-space region over which wavefunction intensity is 
being smoothed. Specifically, theorems by Schnirelman, 
Zelditch, and Cohn de Verdiere (SZCdV) ^ state that for 
a classically defined operator, the expectation value over 
almost all wavefunctions converges to the microcanoni- 
cal average of the classical version of the operator, in the 
7i ^ limit. Since the classical symbol of the operator is 
kept fixed as the limit is taken, these theorems provide 
information only about the coarse-grained structure of 
the eigenstates, and not about the structure at quantum 
nrcchanical scales. 

Wavefunction scarring, the anomalous enhancement 
(or suppression) of intensity near an unstable periodic 
orbit, is a well-known example of non-RMT behavior of 
eigenstates in a classically chaotic system. The distribu- 
tion of wavefunction intensities on a fixed periodic orbit 
can be computed in the semiclassical limit using the lin- 
ear and nonlinear theory of scars &fl|, and is found to 
be very different from the Porter-Thomas prediction of 
RMT. Furthermore, upon ensemble averaging, a power- 
law wavefunction intensity distribution tail is obtained 
(and numerically observed) in chaotic systems, in con- 
trast with the exponential falloff prediction of RMT. The 
fraction of strongly scarred states remains finite in the 
h ^ limit. Nevertheless, scarring poses no threat to 
the SZCdV ergodicity condition, because the size of the 
scarred phase-space region surrounding the orbit scales 
as h, tending to zero in the semiclassical limit. A finite 
intensity enhancement factor affecting an ever smaller 
region of phase space is entirely consistent with ergodic- 



ity on coarse-grained scales. However, the scarring phe- 
nomenon does have very significant effects on physical 
quantities that depend on fine-scale structure, such as 
conductances and decay rates through small (or tunnel- 
ing) leads 0^. 

Another example of markedly non-RMT behavior still 
consistent with SZCdV coarse-grained ergodicity is found 
in the "slow ergodic" systems, such as the tilted wall 
billiard and the sawtooth potential kicked map M. In 
these systems, the classical rate of exploration in mo- 
mentum space is logarithmically slow, and for large 7i" , 
the number of channels occupied by a typical eigenstate 
scales only as h~^''^ log ti~^ , constituting an ever decreas- 
ing fraction of the 0{h~ ) total available number of chan- 
nels. However, the "bright" channels occupied by a given 
wavefunction tend to be evenly distributed over the en- 
tire phase space, and thus coarse-grained ergodicity still 
holds in the limit, even though the wavefunctions are be- 
coming less and less ergodic at the single channel scale 
as ?i — > 0. The present paper extends the indications 
of non-RMT "clumping" of wavefunction density beyond 
the effects of scarring. Moreover, we use as our examples 
the original paradigm of classical Hamiltonian chaos, the 
Sinai billiard, and closely related quantum maps. 

The remainder of this paper is organized as follows: 
in the next section we discuss measures of "microscopic" 
(fi,— scale) quantum ergodicity, including various inverse 
participation ratios and channel-to-channel transport 
measures. Then in Section III the connection is made 
between these stationary properties and the short-time 
dynamics of a quantum system. In Section [V the Sinai 
kicked map, a one-dimensional model for the Sinai bil- 
liard is introduced and discussed. Strong deviations from 
single-channel quantum ergodicity are predicted, and dis- 
tributions for various quantities are computed, that dif- 
fer greatly from RMT expectations. We see that clas- 
sical methods can be used to determine the non-ergodic 
structure of the quantum wavefunctions, even though the 
classical dynamics is entirely ergodic. Quantitative com- 
parison with numerical data follows in Section M. In Sec- 
tion W^ a similar analysis follows for the two-dimensional 
Sinai billiard system, a paradigm of classical and quan- 
tum chaos. Here again strongly non-RMT wavefunction 
intensit y di stributions are predicted and observed. In 
Section VH a simple matrix model is presented and stud- 
ied, some of the statistical properties of which correspond 
to those of the Sinai systems. Similarities and differences 
between the Sinai systems and the matrix ensembles are 
discussed. In the final section we sum up the results and 
discuss certain directions for the future. 



II. MEASURES OF ft-SCALE ERGODICITY 

We now review some important concepts related to 
the quantitative measurement of quantum structure and 
transport at "microscopic" (i.e. single-wavelength or 



single-channel) scales. An alternative discussion may be 
found in j|]. 

Consider a classically ergodic system with quantum 
eigenstates |^) and a test state basis \a). The test basis 
can be chosen to be the set of position states, momen- 
tum states, phase-space Gaussians, or any other set of 
states motivated by the physics of the problem. Often 
the test basis will be taken to be the set of eigenstates of 
a zeroth-order Hamiltonian Hq, of which the full system 
Hamiltonian H = Hq + SH is a perturbation. One is then 
interested in determining whether the true eigenstates |^) 
have a nontrivial structure in the states |a), or whether 
the perturbation SH is sufficiently large so as to com- 
pletely randomize the matrix elements (a|^). Thus, in 
the case of tight-binding models (e.g. Anderson localiza- 
tion), one may consider Hq to be the Hamiltonian with 
on-site energies only, and \a) to be the position states. 
The matrix elements (a|^) then measure the degree of 
localization in position space as the hopping matrix ele- 
ments are turned on. Similarly, in a scattering problem 
one often finds it useful to use momentum states or chan- 
nels as the reference basis |a), and look for localization 
of the full eigenstates relative to this basis. 

For simplicity, we assume that the classical dynamics 
given by H completely mixes the states \a) with each 
other, so that no conservation laws prevent each of the 
eigenstates |i^) from having equal overlaps with all of the 
test states. In the presence of energy conservation or 
other conserved quantities, the formalism outlined below 
needs to be modified to take into account constraints im- 
posed by the classical symmetries. This can be done in 
a straightforward way by, for example, taking the test 
states I a) to be coherent states (Gaussians) in phase 
space. Then it is easy to compute the classical intersec- 
tion of each such Gaussian with any given energy hyper- 
surface, and the actual quantum intensities |(a|^)p can 
be normalized by this classical result. In this way one 
can easily identify the degree of eigenstate localization 
(or deviation from ergodicity) due to quantum effects, as 
opposed to purely classical constraints. See Ref. ||l^ for 
a fuller discussion. 

We will then focus on the set of (properly normalized) 
overlap intensities 
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to devise measures of "microscopic" localization or er- 
godicity in the system under study. In RMT (a natural 
baseline assumption in the absence of dynamical infor- 
mation about our system), the (a|^) are predicted to be 
given by uncorrelated random Gaussian variables, real 
or complex. The intensities Pa^ then follow a x^ dis- 
tribution, of one or two degrees of freedom, respectively. 
Quantum localization will produce an excess of very large 
and very small intensities, compared to this baseline re- 
sult. For convenience, we adopt the normalization where 
the mean intensity is set to unity: 



{Pada = {Pai)i = 1 



(2) 



Here the averages (. . .)(^ are taken over all eigenstates |^): 
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where N is the total number of states accessible from 
I a) (the dimension of the effective Hilbert space). The 
averaging (. . .)a over basis states \a) is defined similarly. 
It is often convenient to compress the intensity infor- 
mation into the set of (local) inverse participation values 
(IPR's) 0: 



IPRa 



{Pl,)^- 



(4) 



Paa is a convenient alternative notation for IPRa, as we 
will see below when we discuss transport in Eqs. M to O. 
For a given test state |a), the IPR at \a) gives the first 
non-trivial moment of the Pas^ distribution (namely, the 
ratio of the mean squared intensity to the square of the 
mean), and thus gives a concise measure of the degree 
of localization at \a) Q The IPR measures the inverse of 
the fraction of eigenstates which have significant inten- 
sity at I a). Thus, equal intensities of all the eigenstates 
at I a) would imply IPRa = 1; this level of ergodicity 
is of course almost never achieved in a chaotic system. 
Gaussian random fluctuations (RMT) produce IPR's of 
3 (for real overlaps (a|<^)) or 2 (for complex overlaps). 
IPR's exceeding the appropriate baseline value signal the 
presence of a localization mechanism beyond RMT. In 
the extreme localization limit where one eigenstate has 
all its intensity at |a), we obtain the maximum possible 
value, IPRa = N . [A prime example of such extreme 
behavior is the case of the "bouncing ball" states [ p^ , 
associated with non-isolated, marginally stable classical 
periodic motion. Such classical trajectories can trap a 
quantum wavepacket \a) for a time comparable to (or 
even longer than) the time at which individual quantum 
states are resolved, causing the wavepacket to have 0(1) 
overlap with only one or a few eigenstates, and leading 
to IPRa = 0(N). From the SZCdV theorems, we easily 
see that that the fraction of bouncing ball states must 
tend to zero in the ?i — * limit. This kind of localization 
is easily visible to the naked eye; other kinds of localiza- 
tion, where the number of eigenstates having intensity at 
the test state \a) is large compared to 1 but small com- 
pared to the total number of states iV, may be less easy 
to detect visually but may also be a statistically more 
important correction to RMT predictions, surviving at 
arbitrarily small values of ?i.] 

One can similarly define an eigenstate-specific IPR in 
the |^)-basis: 



'^A slightly different measure of eigenstate localization at 
a given test state, defined in analogy with classical entropy 
ideas, is discussed in uM- 



where the average is taken over the test basis \a). This 
of course measures the inverse of the fraction of phase 
space occupied by a given eigenstate |^), in the |a)-basis. 
A global IPR can also be defined: 



IPR = (IPRa)a - (IPRc)c 



(6) 



This last quantity measures the inverse fraction of phase 
space occupied by the average eigenstate (or equivalently, 
the inverse fraction of eigenstates that have intensity at 
an average location), and can serve as a simple figure of 
merit for the degree of quantum localization in a given 
system. 

We can relate eigenstate localization to dynamics in 
the following way. Let 



Aaa{t) 



{a\e-'"'\a) 



(7) 



be the return amplitude for state \a) to come back to 
itself after time t. Given a non-degenerate spectrum, the 
mean return probability of the state \a) at long times is 
proportional to its eigenstate IPR: 



Paa = N{\Aaa{t)\'')t 



(8) 



as is easily seen by inserting complete sets of eigenstates 
on the right hand side. Here the time average on the right 
hand side of Eq. ^ is taken over times large compared 
to the Heisenberg time Th, i.e. fi over the mean level 
spacing. 

Similarly, we can relate long-time transport to eigen- 
state correlations. Defining the transport amplitude 



Aab{t) ^ {a\e-'nh) , 
and a cross-correlation analogous to Paa, 

Pab = (PaiPbi)^ , 

we immediately see 

Pab^N{\Aab{t)\^)t, 



(9) 



(10) 



fll) 



where again the time average has been taken of the long- 
time transport probability from \b) to \a). Of course the 
total probability summed over final states for any given 
initial state is normalized: 



(Pab) a = 1 



(12) 



for each b. The simplest non-trivial quantity which will 
measure the fluctuation in the probabilities of being in 
various final states \b) given an initial state \a) (or vice 
versa) is: 



Qa 



{p'D^- 



(13) 



Roughly speaking, Qa measures the inverse fraction of all 
channels that are accessible at long time from channel a. 



In RMT, all the transport probabilities Pab are equal to 
unity with small fluctuations (except for the enhanced 
return probabilities Paa = 2 or 3), so Q^ = 1 for each 
I a) in the iV — > oo semiclassical limit. Qa > 1 indicates 
uneven visiting of the available state space starting in 
the initial state \a), and the overall ergodicity of long- 
time transport can again be summarized in 



Q = (Qa) 



(14) 



In the slow ergodic systems such as the tilted wall bil- 
liard and the sawtooth potential kicked map Q , a highly 
anomalous IPR measure was predicted and observed for 
small h, with the system-averaged IPR (Eq. |g) scaling 
as vh~^/\ogh~^. Semiclassically, the degree of localiza- 
tion in such systems is even stronger, with the IPR scal- 
ing as ?i~^/log?i~^. The difference is caused by diffrac- 
tion, which dominates the phase space exploration and 
increases by vh^ the fraction of phase space occupied by 
a typical eigenstate. These same diffractive effects lead 
to almost perfect long-time transport between channels, 
with the transport measure Q an /i-independent constant 
in the ?i — > limit. 



III. SHORT-TIME DYNAMICS AND THE 
EIGENSTATES 

We now discuss the way in which short-time dynamics 
produces lasting effects on stationary properties, such as 
eigenstate localization and long-time transport |{4|,^,|lO| . 
Define the local density of states (LDOS) at \a) as the 
Fourier transform of the autocorrelation function Aaa (t) : 



N f 
SAE) ^ ^^ j 



dt~ 



iEt 



Aaa{t)^Y.P<^^^-^^^- (15) 



The linearity of the Fourier transform implies that large 
short-time recurrences in Aaa{t) get "burned into" the 
spectrum, producing an envelope which must be the 
smoothed version of the full spectrum Sa{E). Thus, let 



Aaa{t)^Afrit)+A^::^{t) 



long/ 



(16) 



(The most convenient separation between short- and 
long-time recurrences is situation-dependent, as we will 
see below.) Then the full spectrum is given by the sum 
of a short-time envelope and a high-frequency oscillatory 
structure, coming from the long-time dynamics, that is 
superimposed on top of that envelope: 



Sa{E) = ST'\E) + ST^{E) 



(17) 



In the presence of chaos, the number of classical re- 
turning trajectories leading from any state \a) back to 
itself grows exponentially with time (with some Lya- 
punov exponent A). It is then convenient to classify as 
"short time" those returns that are governed by one or 



a small number of classical paths, and as "long time" 
those that arise from interference between many classical 
paths, and for which a statistical description is valid. Due 
to the exponential proliferation, the dividing line between 
these regimes is sharply defined in the semiclassical limit 
7i — !■ 0, being given by the mixing time Tmix ~ j fog ^" ■ 
Of course, no harm would be done were we to err on the 
safe side by computing explicitly intermediate-time re- 
turn amplitudes which could instead have been treated 
statistically. 

Now, in the presence of chaos, the long-time returning 
amplitude at time t is expected to fill the initial test state 
I a) in a uniform, unbiased way |q|, so that the subsequent 
evolution of this newly returned amplitude is equivalent 
to the evolution of the original state. More explicitly, we 
may write for small r 



= (a|a(i))(a(i)|a(M 

= A'::^{t)Af:^'{T) 



r)) 



(18) 



Thus, randomly returning amplitude at long time t leaves 
its imprint on nearby times t + t, and following H we 
may write the full long-time amplitude to return as a 
convolution |j 



A'::^(t) = Y,AZ\t-r)Af°^\r), 



(19) 



where A™*^ is an uncorrelated Gaussian random variable; 



{Al^a 



{t + A)AZ\t))^-5^0: 



(20) 



the averaging being performed over long times t or over 
an appropriate ensemble. The 1/iV factor provides the 
right normalization for the probability to return in the 
absence of nontrivial short-time overlaps, i.e. when 

Af°^\T)^5r,^. ^ . . _ 

Fourier transforming, the convolution in Eq. n^ leads 
to multiplication of the smooth short-time envelope by 
random oscillations in the energy domain: 



S'r>^{E)^Sf°'^'{E)ST''{E). 



(21) 



At very long times (beyond the Heisenberg time Tjf, 
which scales as fi divided by the mean level spacing), the 
spectrum Sa{E) becomes resolved into individual spec- 
tral lines 



^Discrete-time notation is used here for simplicity and be- 
cause it is most useful for the repeated scattering situation dis- 
cussed below. Refer to Ref. H] for a full treatment of continu- 
ous time, which involves introducing an additional very short 
time scale associated with the initial decay of the wavepacket 
I a) (or with Ti divided by the energy uncertainty of \a)). 



Sa{E) 






a^Sr'^\E)SiE ~ E^) 



(22) 



where r^j has the statistical properties of a x^ variable. 
That is, 



rai, 



== 1 



'a^, 



= F, 



(23) 



where averaging may be performed over eigenstates |^), 
test states \a), or over an ensemble, and the constant F 
is given by 2 or 3, for complex or real eigenstates, respec- 
tively. RMT predictions are recovered in the dynamics- 
free limit 5*^'^°'* (£') — 1; short-time recurrences cause 
Sf^™*^{E) to vary with energy. This variation in turn 
leads to larger-than-expected wavefunction intensities 
Pa^ at some energies and smaller intensities at others, 
corresponding to an enhanced IPR and deviation from 
microscopic ergodicity. 

The formalism outlined above has already been used 
to study the statistical properties of the scarring phe- 
nomenon, the anomalous enhancement of certain quan- 
tum eigenstates along the unstable periodic orbits of the 
corresponding classical chaotic system. There, the test 
state I a) is a wavepacket launched on or near the clas- 
sical periodic orbit, and short-time quantum recurrences 
can be computed analytically (for small h) in terms of 
the monodromy matrix and action of the classical or- 
bit. One finds that the IPR for a test state on the orbit 
scales inversely with the instability exponent A of the 
orbit 1^ , and the full distribution of wavefunction inten- 
sities on and off the orbit can be computed as a function 
of A [|6| . Averaging over an ensemble of chaotic systems 
was shown to produce a power-law tail in the intensity 
distribution, dominated by very strongly scarred states, 
and in contrast to the exponential tail prediction of RMT. 
Antiscarred states (ones with anomalously low intensity 
in certain regions of phase space) are of great importance 
in open systems: for example, they have been shown to 
dominate the long-time quantum probability to remain 
in a classically chaotic system coupled to the outside via 
a single-channel lead. The size of the effect is exponen- 
tially large for small A M. More recently, optimal test 
states for measuring scarring ( "universal scarmometers" ) 
have been developed, which take into account an entire 
classical orbit and the linearized classical dynamics in its 
vicinity [|l3| : they provide larger IPR's and more evi- 
dence of wavefunction localization than do simple Gaus- 
sian wavepackets. 



IV. SINAI KICKED MAPS 

A. Definition of system and motivation 

The Sinai billiard uM is a prototypical example of 
strong classical chaos: it consists of a point particle 
bouncing freely in a rectangular cavity with hard walls, 
with a hard disk obstruction placed in the center of the 



rectangle. The system has positive entropy classically 
for any disk size; of course, this fact becomes relevant to 
the quantum mechanics only in the limit where the quan- 
tum wavelength is small compared to the size of the disk. 
(The mixing time Tmix after which a typical wavepacket 
spreads over the entire available phase space is then short 
compared to the Heisenberg time Th, defined as h over 
the mean level spacing, at which the quantum dynam- 
ics becomes quasi-periodic and individual eigenstates and 
eigenvalues begin to be resolved.) 

The statistics of energy levels in the (desymmetrized) 
Sinai billiard has been found to be in good agreement 
with the GOE predictions of random matrix theory [p5| . 
On the other hand, the eigenstate structure of the Sinai 
billiard turns out to be very different from RMT expec- 
tations, and the inclusion of short-time dynamical effects 
is essential for understanding its quantum ergodic prop- 
erties. We will return to a detailed discussion of the 



classical and quantum Sinai billiard in Section VI 



Here we begin with a simplified one-dimensional model 
which contains most of the important structure of the 
original two-dimensional system. We notice first that 
finding eigenstates of a given symmetry class in the Sinai 
billiard is equivalent to finding the eigenstates in a rect- 
angle one-fourth the original size, with a quarter-circular 
bump in one of the corners (and possibly with Neumann 
boundary conditions along one or both of the sides meet- 
ing at that corner). We can then imagine finding the 
eigenstates using an S-matrix approach ||l^,|l^, where 
one considers the scattering of channels of the "free" 
rectangular system off the quarter-circular bump. The 
S-matrix has a strong diagonal component due to the 
straight part of the wall containing the bump, and a com- 
plicated off-diagonal structure due to actual scattering off 
the bump. The long-time dynamics and stationary prop- 
erties of the system (e.g. eigenstates and eigenvalues) are 
obtained by iterating the scattering process. 

From the surface of section method, we know that 
Hamiltonian dynamics in a two-dimensional configura- 
tion space at fixed energy is dimensionally equivalent 
to a discrete-time mapping of a one-dimensional sys- 
tem, and can in fact be reduced to such a system. The 
one-dimensional model we consider in this section is the 
"Sinai kicked map" , defined on a two-dimensional phase 
space {q,p) e [0, 1) X [0, 1): 



p —f p = p — V {q) mod 1 
q —>■ q = q + p mod 1 . 



(24) 



The equations of motion Eq. |2j can be viewed as arising 
from a potential that is periodic in time: 



H{q,p,t) ^ ^ + Viq)J2^it - n) . 



(25) 



At the beginning of every time step, the particle is 
"kicked" by the potential V, following which the po- 
tential is turned off and free evolution takes place for 



a unit time interval. The process is then iterated to ob- 
tain the long-time behavior. The Sinai billiard's straight 
wall with a bump has its analogue in the kick potential 



V{q) 



K 

= otherwise, 



for 



< 



/ 
2 

(26) 



with a parabolic bump (centered at q — 1/2) of spatial 
extent / < 1 . X is a constant (which we will set to be of 
order unity) that determines the typical impulse exerted 
by the bump. The parabolic shape of the potential bump 
is chosen for simplicity only; none of the discussion be- 
low would be affected if a semicircular or other curved 
potential were used instead. The key property of the re- 
pulsive potential V{q) is the parameter /, which sets the 
fraction of an incoming wave that is scattered classically 
after one iteration of the map. 

The Sinai kicked map is a hard chaotic system with 
no stable phase space regions (as can be seen easily by 
computing the Jacobian of the iterated mapping, using 
the fact that V"{q) < everywhere Q). Like the Sinai 
billiard (and the Bunimovich stadium), the system has 
a measure zero set of marginally unstable trajectories, 
given for example hy \q — ^\ > ^, p = 0. After quan- 
tization, such orbits will give rise to "bouncing ball" 
states |12 , which are very strongly localized in momen- 
tum space near p = 0. Our primary interest, however, 
will be not in this measure zero set of states, but rather 
in the structure of the "typical" quantum wavefunctions, 
which obey SZCdV coarse-grained ergodicity, yet have 
very non-uniform structure at the single channel scale. 

The quantization of kicked systems of the form Eq . |2^ 
is straightforward and well-covered in the literature [p_9[. 
A value of h should be chosen so that N = 1/2-kT%, 
the number of Planck cells covering the toroidal clas- 
sical phase space, has an integer value. Then an iV- 
dimensional position basis for the Hilbert space is given 
by qt = [i + to)/N, i = 0...N - 1. Similarly, the 
momentum space basis is given by pj = (j + ei)/iV, 
J = . . . iV — 1. eo,i form a family of possible quanti- 
zation conditions (they correspond to phases associated 
with circling the torus in the p and q directions, respec- 
tively). The two bases are related by a discrete Fourier 
transform. The quantum dynamics is now defined by the 
unitary N x N matrix 



U = exp 



~* 1 Ip' ] /ri 



cxp HV{q)/n] , (27) 



where each factor is evaluated in the appropriate basis, 
and an implicit forward and backward Fourier transform 
has been performed. 

We are now ready to study the structure of the Flo- 
quet or scattering eigenstates of the quantum dynamics 
U, in the pj basis. We notice first that because of the 
symmetry of the kick potential V{q) under q — > —q, the 



classical system has a time reversal symmetry and a par- 
ity symmetry: 



T : t ^ — i, q — > —q 
P : p^ -p,q-> -q. 



(28) 



It will be convenient for us to choose a nonzero value for 
the boundary condition parameter ei , thus breaking the 
parity symmetry P under quantization, while maintain- 
ing the time reversal symmetry T by setting ep = 0. The 
eigenstates are then real in the momentum basis, and 
the appropriate RMT baseline is IPR = 3 (see discussion 
following Eq. 0). For an asymmetric bump or kick po- 
tential, the quantum wavefunction intensity fluctuations 
would be expected to follow a x^ distribution of two de- 
grees of freedom under RMT, giving rise to the baseline 
value IPR = 2. The analysis to follow is of course com- 
pletely independent of the symmetry chosen, provided 
that the appropriate baseline quantum fluctuation factor 
F is used. 



B. Short-time dynamics 

1. Quantum factor of two 

As suggested in the preceding section, we should be- 
gin our analysis by examining the classical and quantum 
short time dynamics of the Sinai kicked map in momen- 
tum space. In analogy with the quantum return ampli- 
tude Ann{t) of Eq. 0, let Pnn{t) bc the classical probabil- 
ity to remain in state p„ after t iterations of the map. ^ 
Classically, for any incoming momentum p„, a fraction 



Plil) = l-f 



(29) 



of all particles remain in momentum p„ after one kick, 
while the remaining fraction / get scattered to other mo- 
mentum states. Semiclassically, there is an amplitude 
1 — / for remaining in the incoming channel, as can eas- 
ily be seen by taking the (semiclassical) overlap of the 
initial and final states. The probability to remain unscat- 
tered after one step is then 



\Unn\' = \Annil)\' = {I - ff 



(30) 



Notice that the quantum one step survival iC/nnp is 
smaller than the classical probability P^ni^) f^^' '^ot scat- 
tering. In analogy to ordinary scattering in free space, 
one can define a cross section for scattering off the de- 
fect. In the limit of small /, we see from the above anal- 
ysis that the quantum cross section is twice as large as 



^For conventional resons we will be using the index n to label 
the momentum states instead of the generic index a used in 



Section III 



the classical. (This is in complete analogy with quan- 
tum scattering theory in free space, in which diffrac- 
tion results in a quantum cross section twice as big as 
the classical, even in the short wavelength limit. Essen- 
tially, the far-field diffraction into the shadow zone dou- 
bles the quantum cross section). Here, there is a quantum 
probability /(I — /) for diffracting into nearby channels. 
Classical-quantum correspondence still holds after appro- 
priate coarse-graining over scales large compared to h/ f 
(but still small classically) in momentum space. As men- 
tioned earlier in our discussion of the Sinai billiard, we 
are always working in the semiclassical regime fi <^ f , 
where the bump size is large compared to a wavelength, 
though it may be small compared to the system size 1. 

This difference between the classical and quantum 
probability to remain in the initial channel ((1 — 
/) vs. (1 — f)"^) survives the limit ?i — > 0. In this limit 
the fraction of diffracted amplitude and the number of 
scattering channels into which diffraction occurs both re- 
main constant. For finite ?i, there will of course be an 
O^yK) correction to the quantum amplitude A„„(l), as 
it is possible for the horizontal {V'{q) = 0) portion of the 
bump to scatter an incoming channel p„ back into itself. 



nearby channel in one step and diffracts back into p„ 
during a following scattering event. The extra amplitude 
coming from diffracting back and forth between nearby 
channels can add in or out of phase with the "naive" 
contribution. As we found previously, the probability to 
scatter into a nearby channel after one step is /(I — /); 
this is comparable to the probability / for scattering into 
a classically distant channel and being completely lost 
from the system as far as the short-time return probabil- 
ity is concerned. 

Notice that the short-time return probability of Eq. pll 
is completely independent of the shape of the bump V[q) 
(as long as the bump amphtude K is chosen to be 0(1) 
so as to allow scattering into many distant channels). 
In fact, we can compute the short time return proba- 
bility in a simplified model where the non-zero part of 
the potential is replaced by an absorber, and pieces of 
the probability density simply get subtracted from the 
system. For given i, we then have a distribution of the 
quantum probabilities to remain |A„„(i)p. We easily see 
that the fastest possible decay of the initial state p„ is 
obtained for p„ = /, where an entirely untouched piece 
of the wavefunction is absorbed at each step: 



2. Multistep scattering 

We proceed to analyze the multi-step behavior of the 
dynamics, particularly the probability to remain in the 
initial state p„. In the absence of step-to-step correla- 
tion, the classical probability to remain unscattered af- 
ter 2 steps would be P,'jjj"^'^<=(2) = (1 - /)^ giving rise 
to a quantum probability |A;^?f™(2)|2 = (1 - /)■*. This 
is the same probability that we would obtain by simply 
iterating the diagonal part of the evolution matrix, i.e. 

by approximating (f/^)„„ = J^n'^nn'Un'n « Un„Unn- 

Of course, the true probability to remain in state p after 
two steps is p-dependent: for most values oi Pn, namely 
\Pn\ > /, entirely different parts of the wavefunction are 
scattered at each of the two steps, so the classical prob- 



abihty to remain is P,f„{2) = 1 - 2/ < P; 



cl, naive 

nn 



(2). On 



the other hand, for p very close to zero, most of the 
probability that would be scattered at the second step 
has already been lost on the first scattering event, so 
P„=U2)«l-/>P,^k"^'™(2). 

The analysis can be extended easily to longer times. 
The quantum probability to remain after t steps is still 
given, to leading order order in h, by 



\{U%n\'^\AUt)\' = \Prfn{t)\\ 

where naively (in the absence of correlations) 



^cl, naive 



{t) = {i-fY. 



(31) 



(32) 



The true value of P^^(i) for i > 1 is p-dependent; in 
quantum mechanics, this p„— dependence can be under- 
stood in terms of amplitude that diffracts from p„ to a 



p:^(t) = i-tf^ 



(33) 



for t < 1//. The largest values of P^l^ arise from p„ near 
zero, as described above, and also from p„ that are near 
simple fractions like 1/2 or 2/3. These slowly-decaying 
momentum channels give rise to the most non-ergodic 
long-time quantum behavior, as we shall see below. 

For each channel p„ , then, we can compute the quan- 
tum short-time autocorrelation Ann (t) ; it is given by the 
square root of the quantum probability to stay (Eq. ^l|) , 
times the phase accumulated from the (free) quantum 
dynamics: 

Atnt) = iPnmPn) = e-''"*/2p„t (1^1) • (34) 

This holds for both positive and negative short times 
(note A{—t) — A*{t) by unitarity). For a typical mo- 
mentum Pn, A^™*'{t) has a decay time of 0(1//); upon 
Fourier transforming we obtain a short-time spectral en- 
velope S^^°'^*'{E) centered at quasienergy p^/2 and with 
width of order /. Specifically, using the naive estimate of 
Eq. ^ and taking the bump size / to be small, we obtain 
a Lorentzian short-time envelope 



for 



SI 



E-E^ 



short,] 



'{E) = 



2/ 



P 



P-^ 



(35) 



<1. 



C. Long-time behavior and stationary properties 



1. Scaling properties 

At times long compared with 1//, most of the initial 
amplitude in a typical channel p„ will have been scat- 
tered by the bump, and the return amplitude A^°^^{t) 
will be given semiclassically by a sum over many nontriv- 
ial paths (the relative phases between the paths being of 
course ?i— dependent). As discussed in Section III, and in 



completely analogy with nonlinear scarring, these long- 
time recurrences are given by independent Gaussian ran- 
dom variables, convoluted with the short-time dynamics 
Af°'\t). The full local spectrum Sn{E) (Eq. |l5|) is a 
line spectrum with individual intensities P„j = A^K^^IOP 
(where we have denoted {p„\ as (n|) being given by a x^ 
distribution, weighted by the height of the linear envelope 
at energies E^: 



M 



■a^ST'\E^) , 



(36) 



Here r^^ are independent x^ variables with mean unity 
(see Eq. HI) . The expected local IPR (Eq. ^ is given by 
a product of a factor associated with the short-time en- 
velope and a factor {F = 3) associated with the spectral 
fluctuations r^j under the envelope: 



-t-oo 



IPR„ = 3((5f -*(i?)) )b - 3 ^ |At"'(OP • (37) 



Using the naive short-time dynamics of Eq. 
an estimate for the typical IPR: 



IPR" 



3x 



2 - 2/ - /^ 

V-P 
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we obtain 



(38) 



Using the upper bound we obtained in Eq. |3^ on the rate 
of short-time decay of an initial momentum channel, we 
also have a lower bound on the local IPR: 



IPR" 



"^3/-/ 



(39) 



Notice that this lower bound is for moderate / already 
larger than the RMT expectation of 3. 

For a given value of /, we may use our knowledge of the 
short-time classical dynamics and Eqs. |4[ ^ to obtain a 
distribution of the local inverse participation ratios IPR„. 
For / ^ 1, the decay time of the classical autocorrelation 
function P^J^(i) scales with 1//, so we expect the IPR 
distribution Vj to scale likewise: 



7'/(IPR„ = x) = fVifa 

for some function V{y). From Eq. t 
bound on possible IPR's for small /, i.e 

r{y) =Ofory < 2. 



(40) 
we have a lower 

(41) 



Using the naive estimate of Eq. pq for the IPR at a "typ- 
ical" value of the momentum, we determine that the me- 
dian of the distribution Vly) should be in the neighbor- 
hood of 3. This discussion of the IPR distribution has 



been very general; however the details of the function 
P{y) may in fact depend on classical system parameters 
other than the bump size /. For example, in the equa- 
tions of motion Eq. g^ we could have replaced the free 
evolution in the second line with 



q 



q = q + ap mod 1 



(42) 



making the elapsed time a between kicks an arbitrary 
parameter. (In the Sinai billiard system, the parameter 
a corresponds roughly to the length-to- width ratio of the 
rectangular billiard.) The detailed properties of the IPR 
distribution V{y) will then depend on the values of clas- 
sical parameters such as a, while results such as Eq. ^ 
apply more generally to the entire class of Sinai-type sys- 
tems. Below, in Figs. U and 0, we present the actual 
classically computed function V{y) for the Sinai kicked 
map with a — 1; there Viy) is compared with statistics 
collected for the corresponding quantum system. 



2. Tail of the IPR distribution 

First, we discuss another important qualitative feature 
of the IPR distribution, namely the long tail oiP{y) com- 
ing from momentum channels Pn which are near simple 
fractions and thus decay on a time scale longer than the 
typical 0(1//) steps. Consider a very small initial mo- 
mentum \p„\ <^ f. As we saw in the discussion imme- 
diately preceding Eq. |3l|, only a very small fraction of 
the remaining amplitude in \pn) is scattered during each 
kick following the first one, because the part of the wave 
which has not yet been scattered shifts very little in posi- 
tion space between kicks. Explicitly, the classical proba- 
bility (and thus the quantum amplitude) to remain after 
t steps is given by 



i^?r(oi=^i(i) 



l-f-{t- l)Pr, 



(43) 



for 1 < i < (1 — f +Pn)IPn- Thus the decay time for the 
initial state \p) scales as 1/pm and the inverse participa- 
tion ratio IPR„ scales likewise (compare Eq. B^. These 
IPR's are large compared to the 0(1//) IPR's obtained 
for typical channels (yet small compared to the 0{N) 
IPR's which characterize bouncing ball states). 

Similarly, if we choose a momentum channel which lies 
near a simple fraction, p^ — ^ + e (for le < f < l/£), 
then following the first £ kicks, a fraction £e is scattered 
after each successive kick, and the decay time (and IPR) 
for such a channel therefore scales as l/£e. We dub 
these special momentum channels the "near-bouncing 
ball" trajectories. We can now easily estimate the frac- 
tion of channels with IPR greater than some number x, 
where x ^ 1//. All channels within 1/x of zero sat- 
isfy this condition, as do those within l/£x of a simple 
fraction m/£. Now we note that for a typical integer 
£, a finite fraction of integers m = I . . . £ are relatively 
prime to £; thus from each value of £ we obtain a fraction 



0(£ X j^) = 0{l/x) of channels with IPR expected to be 
greater than x. Adding up contributions from all values 
oil between 1 and 1//, we have a cumulative probability 
0(1/ fx) for IPR„ to be greater than x, or 



P(IPR„ = x) - ^ 



1 
/~2 



(44) 



for X ^ 1//. We see that the parameter / enters in 
the expected way, and the tail of the scaling distribution 
(compare Eq. ^^ is then given by 



for y :$> 1. 



V{v) ^ l/y2 



3. Failure of channel level ergodtcity 



(45) 



From the inverse square form of the IPR distribution 
tail in Eq. |4J, it would appear that the mean value of 
the IPR diverges logarithmically for these systems. How- 
ever, we notice that at fixed energy the possible IPR is 
bounded above by the total number of channels N (this 
being the IPR for a pure bouncing ball state), and so we 
have 



(IPR„)„ = {Pnr. 



logN 



(46) 



We see from Eq. E^ that the mean inverse participation 
ratio in the kicked Sinai systems diverges logarithmically 
with increasing energy (or decreasing h); thus the wave- 
functions are becoming less and less ergodic at the single- 
channel scale even as the classical limit is approached, 
despite the ergodicity of the corresponding classical me- 
chanics. The situation is more surprising here than in the 
slow ergodic systems 0, as in the present case the Lya- 
punov exponent is positive and the Sinai billiards have 
long been considered a prototypical example of strong 
classical ergodicity and mixing. We also note that the 
logarithmically increasing mean IPR in Eq. Eq is due 
not to the bouncing ball states (the fraction of these 
scales as 1/N and thus their contribution to the mean is 
iV- independent) , but rather to the "near-bouncing ball" 
channels, whose decay time is large compared to the 
typical decay time 1// but still small compared to the 
Heisenberg time N at which individual eigenstates are 
resolved. Each such channel contributes to many eigen- 
states of the system, but only a small fraction of all the 
available eigenstates. 

Having made predictions about the structure of the 
IPR distribution for Sinai-type systems (2// lower cut- 
off, 0(1//) median, 0( °^, ) mean, inverse square tail), 
we now proceed to perform a similar analysis for the 
other statistical quantities discussed in Section O. As 
discussed previously, the details (factors of order one) of 
the various distributions and statistical averages will be 
system-dependent, and can be computed explicitly for 



any specific Sinai-type system (as we will do in the fol- 
lowing section). What we are interested in here is the 
universal scaling behavior of wavefunction structure with 
the bump size / and the wavelength 1/A^. 

We consider first the individual wavefunction intensi- 
ties Pn^ in the momentum basis. For the typical mo- 
mentum pm we have seen that the smoothed spectrum 
5,^^™* (-E) has a Lorentzian peak of height scaling as 1// 
and width scaling as f, centered on the optimal energy 
En = ^Pn (see Eq. ^5|). Far from i?„, the smoothed 
spectrum levels off to a height of order /, leading to the 
typical behavior 



pmedian 



/• 



(47) 



Notice that the median intensity is much smaller than 
the mean (cf. Eq. 50). 

Even for the most anomalously localized channels, the 
minimum value of S!^°^'^{E) never falls below 0{f); the 
smallest values of P„^ must therefore arise from x^ fluc- 
tuations multiplying this typical intensity. For complex 
wavefunction amplitudes (p|^), this implies 

P(P< = ^) ^ } cxp(-x//) [x « /] , (48) 

and a corresponding expression is obtained in the real 
case, where P„^ is a x^ variable of mean 0{f) and one 
degree of freedom: 



V{Pni ^X) [X < /] 

^y2^TXJ 



(49) 



The mean of the intensity distribution is of course fixed 
by normalization: 



(Pnf) = 1 



(50) 



where once again the averaging (• • •) can be thought of 
as an average over eigenstates |^), momentum channels 
\pn)i or over some ensemble of Sinai- type systems (where 
e.g. the shape of the bump can be varied while preserving 
its total size /). 



4-. Tails of intensity and transport measures 

Lastly, we turn to the tail of the intensity distribu- 
tion, which we expect to result from large values of the 
smoothed spectrum S^^°'^*'{E). (Fluctuations of the full 
spectrum S{E) around its smoothed value have a x^ 
form. The probability of obtaining an intensity P„^ large 
compared to the short-time prediction S'^^°''*(£'„) is ex- 
ponentially small.) As previously discussed, a fraction 
0{1/ fx) of all momentum channels \p) have a peak in the 
spectrum S^°'^^{E) of height > x, and the width of such 
a peak is then 0{l/x). Therefore a fraction 0{\/ fx"^) of 
all intensities Pn^ satisfy the condition 5'^^°''*(i?„) > x, 
and since the fluctuations in P„j around this smoothed 
value are of order unity, we obtain 






1. 



(51) 



The eigenstate-basis IPR measure IPRj = P^^ (Eq. 0), 
which measures the inverse fraction of channels in which 
a given eigenstate lives may be studied in a manner very 
similar to the channel-basis measure IPR„ = Pnn- From 
Eq. ^, we already know the mean value of the IPR^ 
distribution: 



(IPR?) 



(%) 



log TV 



(52) 



We proceed to study the structure of the "typical" wave- 
function |i^). From Eq. p\ we know that given some 
eigenstate |^), the probability that is has intensity > x 
at any particular momentum \p) is 0{1/ fx^). If we as- 
sume the overlaps of |^) with the different momentum 
states to be uncorrelated, and notice that there are a to- 
tal of A^ momentum channels to overlap with, we see that 
for X < y'N/f there will generically be at least one mo- 
mentum \pn) such that P„^ > X. We now compute the 
contribution to P^j from all intensities P„^ = x between 
1// and ^/N/J: 



P 



«? 



> 



{Pli) 



1// 



dx x^ P(P„e = x) 



dx x 



fx^ 



logN 



(53) 



(recalling N ^ I//)- Thus we see that not only the 
mean, but also the inverse participation ratio for the typ- 
ical wavefunction tends to infinity in the classical limit: 



ipr; 



medi; 



logA^ 



(54) 



The tail of the IPR„ distribution arises from the rare 
intensities P„5 » ^JN/ f, and using Eq. [5^ is easily seen 
to take the form 



^(IPRj ^x)^—^ 



1 



(55) 



Having analyzed the statistical structure of individual 
wavefunctions in Sinai- type systems, we can now proceed 
to examine the quantum transport behavior. Specifically, 
we focus on the long-time transport probability P„„' be- 
tween two channels p„ and p„' , as introduced previously 
in Eqs. |l^, |l^. For two typical channels Pn and p„', each 
of the two smoothed local densities of states has the form 
of a peak of height 1// and width / centered around some 
energy, and then falls off to a value of 0{f) far from that 
energy (see Eq. p3). Since the two peaks are generically 
centered at different energies, \En — En'\ ^ /, we easily 
see that the overlap between the two envelopes is 0(/): 



^median 



/■ 



(56) 



Of course we also know the mean value of this distribu- 
tion by construction: 



(P,, 



1 



(57) 



Large values of the transport measure Pnn' arise from 
those Pn and pn> for which the two spectral envelopes 
S^°'\E) and S^^°'\E) are both anomalously taU and 
narrow, and which also have significant overlap with each 
other. Explicitly, in order to obtain a value P„„' > x 
for large x, we require IPR„ > x, IPR„' > x, and also 
\En — En'\ < -■ The combined probability for these three 
unlikely events scales as 



so 



V{Pnn' > X) 



V{Pn 



1 



1 



1 



fx"^ fx^ 



px^ 



kc> 



f' 



(58) 



(59) 



This a very quickly decaying tail compared to the one ob- 
tained previously for the inverse participation ratio P„„ 
(compare Eq. ^J); thus transport efficiency for this class 
of systems is much less anomalous than the structure of 
individual wavefunctions. This makes sense intuitively 
and is also consistent with the findings for slow ergodic 
systems in Ref. |g]. 

Finally, the remaining measure we must consider is 
the final-state-averaged transport efficiency Qn for initial 
state Pn (introduced in Eq. |l^; see also Eqs. |lO[ |ll| for 
the definition of Pnn' , the long-time probability of getting 
to channel pn' from channel p„). The quantity (Qn), as 
well as the typical value of Qm will be dominated by the 
Lorentzian envelopes governing typical intensities P„j: 



1 ^ 

— -TT 2_^ Pne,Pn'(, 



2/ 



2/ 



dE 

2^ \p + {e~pII2y) \p + {e-pI,I2Y 

4/ 



Ap + {pl/2-pl,l2y 



(60) 



where in the second line we have inserted the typical 
intensity in channel p„ of a state |^) with energy E (from 
Eq. H). Now 



" ~ TV ■^ ^nn' 



d{plm 



4,/ 



2^ \Ap + {pl/2~pl,/2Y 
The mean and the median both scale as 



1 
2/ 



(61) 



{Qn) - QT. 



median 



1 

T 



(62) 



and furthermore in the classical limit A^ — > oo it is ex- 
ceedingly difficult to obtain values of Qn either small or 
large compared to 0(1/ f). For almost any initial chan- 
nel Pm the fraction of final channels p'n to which one can 
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FIG. 1. The distribution of the inverse participation ratio 
IPR„ = Pnn is plotted for bump size / = 0.1. From left 
to right, the four solid curves represent data for A'' — 100, 
200, 400, and 1600. The theoretical prediction (see text) is 
represented by dashes, while the random matrix prediction is 
IPR = 3. 



be transported at long times is 0{f) <^ 1. Bouncing- 
ball (free propagation) channels of course have even less 
coupling to other momentum states (roughly speaking, 
they couple to themselves only, P„„ = 0{N) and thus 
Qn = 0{N)), but these constitute a vanishing fraction 
of all channels in the classical limit. 



V. NUMERICAL TESTS IN SINAI KICKED 

MAPS 

We proceed to a numerical study of the structure of 
wavefunctions in the Sinai kicked systems, focusing on 
those statistical properties which we have treated the- 
oretically in the preceding section. We begin by con- 
sidering the distribution of inverse participation ratios 
IPR„ (Eq. H), each of which measures the inverse frac- 
tion of eigenstates having significant intensity at some 
momentum channel Pn- The bump size / is fixed at the 
moderate value of 0.1, which is small compared to the 
system size of unity, yet large compared to wavelengths 
1/A^ < 0.01 which we are going to consider. In Fig. ||, 
the IPR distribution 7'/(IPR„) is plotted (solid curves) 
for several values of the quantum wavelength: N = 100, 
200, 400, and 1600. In each case, an ensemble has been 
constructed by varying the bump height parameter K 
in Eq. g6|; each realization with K — 0(1) is expected 
to have the same wavefunction statistical properties, as 
discussed in the preceding section. 

Also plotted as a dashed curve in Fig. n^ is a classi- 
cal prediction for P/(IPR„). This quantity is obtained 
by taking a random sample of initial momenta p„, and 
for each of them computing classically the probability 
Pnn{t) to remain unscattered after t bounces. Given the 
short-time classical behavior P^J^(i), we use Eqs. |4[ ^ 
to predict the expected quantum IPR for that momen- 
tum channel p„, eventually leading to the distribution 
shown by the dashed curve. Of course, this is a semi- 



classical {N -^ oo) prediction; in particular, it ignores 
fluctuations in the IPR which result from summing over 
a finite number of eigenstates in Eq. ^ (even in RMT, 
fluctuations around the mean value of 3 are expected for 
finite N). 

Indeed, we see in Fig. |^ that the quantum IPR distri- 
bution does appear to approach the classically predicted 
distribution as N gets large; the convergence with N will 
be studied more quantitatively below in Fig. 1^. By the 
time we reach TV = 1600, the peak of the distribution 
is within 10% of the classically expected value, and is 
shifted by a factor of seven from the naive random ma- 
trix prediction. We also notice that all the IPR's in our 
sample are larger than the value of 3 predicted by ran- 
dom matrix theory, and most are larger by a factor of at 
least five: this is unmistakable evidence of strong devia- 
tions from microscopic quantum ergodicity in the kicked 
Sinai systems. 

Next, in Fig. g we fix the total number of channels at 
N = 1000, and study the IPR distribution for various 
values of the bump size /. Guided by the predicted scal- 
ing relation of Eq. Eo, we choose to plot the distribution 
of the scaled quantity / • IPRn for each value of bump 
size /. For each of / = 0.1 (dashed curves) and / = 0.2 
(dotted curves), two distributions are plotted: one for 
the original kick potential of Eq. g6|, and the other for a 
modified kick potential 



Vi,) = -^ 
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(63) 



The latter potential has a kink at q = 1/2, causing a dis- 
continutiy in the classical dynamics. We see from Fig. 
that the choice of kick potential (Eq. p^ or Eq. ^) has 
no significant effect on the IPR distribution, as long as 
the bump size / is fixed, confirming the universality pre- 
dicted in the previous section. In particular, we notice 
that the flat part of the potential [V'{q = 1/2) = 0] in 
Eq. pq, which scatters any incoming channel back into 
itself, has no discernible effect on quantum localization 
at the energies under consideration. 

The classical prediction for the scaling distribution 
V{f ■ IPR) is also plotted in Fig. for comparison (see 
solid curve). We see very good agreement among the 
4 sets of quantum data at f — 0.1 and / — 0.2; sim- 
ilar scaling behavior with / is observed for the billiard 
system in Fig. ^. Again, the slight discrepancy (around 
10%) between the numerical data and the classical pre- 
diction may be attributed to the finiteness of the energy. 
At these energies, the minimum observed value of the 
IPR appears to be near 1.5//, in contrast to the 2// 
semiclassical limit prediction of Eq. p3. 

In the tail, we predict (Eq. ^^ the inverse square be- 
havior V{x) « c/ fx^ for the IPR distribution, where the 
constant c can be determined to be 0.6 through a detailed 
classical analysis of this system as described above. The 
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FIG. 2. The distribution of the scaled inverse participation 
ratio /-IPR is plotted for / — 0.1 (dashed curves) and / = 0.2 
(dotted curves). A'^ = 1000, and for each value of / two dis- 
tributions are plotted corresponding to different bump shapes 
(see text). 
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FIG. 3. The tail of the IPR distribution is plotted for 
N = 1000, with / = 0.1 (solid curve) and / = 0.2 (dashed 
curve). The theoretically predicted 0.6/a;^ behavior (Eq. HH) 
appears a dotted line. 

c/y"^ tail for y = /-IPR is indeed observed in Fig. |^, where 
the prediction appears as a dotted line on the log-log plot, 
vifhile the solid and dashed curves represent / = 0.1 and 
/ — 0.2, respectively. This data was again taken for 
N = 1000, and we see the power-law behavior persist to 
IPR's of about 300, where the IPR becomes comparable 
to the total number of channels and the theory naturally 
breaks down. 

We recall that this breakdown of the inverse-square law 
at IPR ^ N leads to the prediction of mean IPR growing 
logarithmically with N: (IPR) ~ 0.6 log A^ + const (see 
Eq. [46| ). This behavior is indeed observed for N ranging 
from 100 to 1000; we omit the figure here because an 
analogous plot for the billiard system appears in Fig. |l0| 
in the following section. The median Pnn shows no such 
increase with N; it saturates at « 2.35// independent 
of N (see Eq. p8i). The median IPR for an eigenstate 
(P^j), on the other hand, does grow logarithmically with 
iV, but only half as fast as the mean, in agreement with 
Eq. ||. 

We next turn to the distribution of individual wave- 
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FIG. 4. The distribution of intensities P„{ for / — 0.1 and 
N — 1000 (solid curve). The classical prediction (see text) 
follows the data very closely (dashed curve). Asymptotics 
form for the head (Eq. ^) and tail (Eq. |5lj) of the distri- 
bution are both drawn using small dashes. For comparison, 
the Porter-Thomas distribution of RMT appears as a dotted 
curve. 

function intensities. In Fig. H the distribution of intensi- 
ties P„^ is plotted for / = 0.1 and A^ = 1000 (solid curve). 
The classical prediction (obtained as described in the dis- 
cussion of Fig. 1^ above) is plotted as a dashed curve; the 
difference between data and prediction is barely visible 
except in the very tail where the statistical uncertainty 
in the data becomes important. The two analytic asymp- 
totic expressions: exp{—x/2f)/y/2Trxf for small intensi- 
ties X (Eq. Boh and \/ fx^ for large x (Eq. pw are also 
shown in Fig. IJ. These two expressions are valid for 
X <C / and X :s> 1//, respectively. By contrast, the RMT 
Porter-Thomas prediction (dotted curve) does not agree 
with the data in the head, body, or tail of the distribu- 
tion. See also Figs. ^ and |l^, which focus separately on 
the head and tail of the intensity distribution for the bil- 
liard system, and again find good agreement with theory 
and disagreement with RMT. 

The distribution of transport measures Pnn' has also 
been studied and observed to possess a l//^a;^ behavior 
for X :$> 1//, as predicted in Eq. ^. This data is omitted 
here as very similar behavior is obtained for the billiard 
in Fig. O. The overall transport efficiency Q has also 
been studied and follows the predicted scaling Q ~ 1// 
of Eq. p3, so that only a fraction 0{f) of all channels are 
quantum mechanically accessible at long times starting 
in any one initial channel. 



VI. LOCALIZATION IN SINAI BILLIARDS 

The Sinai billiard was the first nontrivial dynamical 
system shown to be ergodic with positive Lyapunov ex- 
ponent pQl- In this sense it is the paradigm of chaos. It 
is also a unit cell of the Lorenz gas, a periodic array of 
hard disk scatterers (see Fig. ||a). 

For numerical reasons we investigate a modified Sinai 
system with the circular disk off center and jutting only 
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FIG. 5. a) The Lorenz gas and two choices for a funda- 
mental domain, b) A Sinai hke billiard related to the Lorenz 
gas 
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FIG. 6. The modified Sinai system, with a partial disk 
occupying a variable fraction / of the right hand vertical wall. 



part way into the billiard; this is sill a chaotic system 
(see Fig. |b). 



A. Scattering method 

A scattering system closely connected with both the 
Lorenz gas and the Sinai billiard puts the Sinai disk at 
the end of a corridor of length a (Fig. ^) . The scattering 
wavefunction can then be expanded as 



^{x,y) 




iknX 



sm{n7Ty/h) 



^n„'e-2*'=""^e^""'^ sin(nV2//6) (64) 



where for later convenience we have factored out a phase 
exp(2iA:„'a) from the n'-th column of the S-matrix. (If 
there is no scatterer on the right hand wall, this makes S 
the diagonal unit matrix, assuming Dirichlet conditions 
there). Now suppose that we reflect the scattered wave 
from the left wall back towards the right hand side, in 
accordance with the closed billiard problem we wish to 
solve. This can be done by imposing a boundary condi- 
tion at the left wall, which need not necessarily be Dirich- 
let. (We indicate this by using a dashed line to represent 
this wall in Fig. o.) If the wave is reflected from the left 
wall at a; = 0, it returns with a new phase exp(j0) given 
by the boundary condition at the left wall. We define 



Unn' = Snn' exp{-2ikn'a + 14 

Setting ■(/;„ = exp(— ifc„a;) sin(r^7^J//&)/^ 
coming (right-moving) wave is then 



(65) 



fc„, the net in- 



(l + C/ + C/2 + ...)V'n = 



1 



\-u 



Ipn 



(66) 



(see Fig. I). 

Evidently, a bound state can be built up in the bil- 
liard if U has an eigenvalue -1-1. We can diagonalize the 
LZ-matrix and consider the properties of its eigenstates. 
Since [/ is a unitary matrix, its eigenvalues lie on the 
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unit circle. As we change the phase shift (j) ^^ the left 
wall, the eigenvalues will correspondingly rotate around 
the unit circle; each of the N eigenvalues of U (assuming 
there are N open channels) will pass through +1 for some 
(/), so that every eigenstate of U is an eigenstate of the 
closed billiard with some boundary condition at the lefi 
wall and Dirichlet boundary conditions elsewhere. 

If one is willing to search through ranges of energies oi 
of box lengths a one can find a set of eigenstates satisfying 
a particular boundary condition; this is a way of finding 
eigenvalues and eigenstates of the billiard with Dirichlet 
boundary conditions; they are given by eigenstates of U 
with eigenvalue 1 [0,pll. However here we do not seek 
the Dirichlet solutions, since they are not special as fai 
as their localization in the channel space (this has been 
tested numerically). This is of great value in gathering 
the statistics needed here. 

Two typical eigenstates of the [/-matrix are shown in 
Fig. 0; these show fairly obvious non-statistical mixing 
of different directions of propagation in the billiard (non- 
mixing of channels in the scattering approach). 



1. Patterns in the channel (momentum) transport 

A new twist arises in the channel transport measure 
Pnn' 1 which we now present. Heretofore we have been 
making the point that Gaussian random wavefunction 
statistics are much stronger than required for SZCdV er- 
godicity, and that much coarser randomness can still lead 
to ergodic transport classically. Below, we see that trans- 
port in momentum space may even be highly organized, 
but in a way that still permits coarse grained SZCdV 
ergodicity. 

The density plot of the transport measure Pnn' for a 
typical case (/ — 0.1 with 280 open channels, side lengths 
equal) appears at the top of Fig. R A pronounced fringe 
pattern is evident. This pattern changes with the length 
of the billiard, and as we now show represents alternating 
constructive and destructive interference due to the phase 
factors exp[2iA:„a] where a is the box length and A:„ is the 
horizontal wavevector. The S'-matrix itself shows none oi 
this fringing, but it is strongly evident already for S^: We 
have 



nn' — / ^ ^nn" ^n"n' : 



(67) 



and since S is diagonally dominated, the major contri- 
bution to Sf^j^, for n y^ n' is 



^nn'Jnn' ' ^nn' ^n' n' 



i>nn'(i>nn ~^ ^n'n') ("°) 



Of course Snn and Sn'n' can interfere; these diagonal 
elements have factors exp[2i/c„a] and exp[2zA:„'a], respec- 
tively. Subsequent iterations reinforce this interference 
and give very sharp preferred channels that one can end 
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FIG. 7. Two typical eigenstates of the S-matrix for the 
Sinai-like scattering system. 
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FIG. 8. Top: long-time transport probability between 
channels n and n' for the a = b = 2n (square) billiard, f=0.1 
(m = h = 1). Bottom: the fringe pattern from Eq. g8|. 



up in when starting from a given initial channel. A plot 
of 

Wnn' = \ex.p[2ikna] + exp[2ifc„/a]| 

appears at the bottom of Fig. g, and is seen to bear a 
close resemblance to the fringe pattern in P„„/ (the ex- 
ponent 28 is of course arbitrary and only serves to set 
the contrast ratio of the plot). It should be kept in mind 
that the fine detail (pixel by pixel) of the intensity mod- 
ulations present in P„„' are absent in the lower plot, but 
the overall modulation of the regions of large and small 
Pnn' are almost identical. 

Interestingly, the special channels which correspond to 
classical free motion (never hitting the obstruction) show 
up on the diagonal as hyperbolic points of high density. 
This may be shown by expanding in the channel index 
(at least in the lower n region where the Taylor series 
holds for An ^ 1), e.g. 



|exp[2zA:„a] -|- exp[2ifc(„+An)a] 
|l-|-exp[2Ja (dkn/dn) An]\ , 



(69) 



where a is the length of the rectangular box. Since kn — 



'2{E— "2^2 ) and kny = mr/b, where b is the 
height of the box, we have 

dn k„ b 



(70) 



Then the interference in Eq. p9^ is maximally constructive 
for 



k 

'^ny 
kn 



m b 
An a 



(71) 



i.e. exactly for the free motion trajectories. The special 
channels correspond with the hyperbolic regions along 
the diagonal in Fig. @. The near-bouncing ball channels 
near the free propagation channels preferentially diffract 
symmetrically about these special channels, as evidenced 
by the local hyperbolic structure. This is again a con- 
sequence of the interference structure in Eq. ^. Essen- 
tially, there is a preference to scatter by a multiple of a 
reciprocal "lattice" vector, (2aAk = 2rmr), reminiscent 
of Bragg scattering from a periodic structure with lattice 
constant a. 

The dramatic interference pattern is another interest- 
ing quantum signature of a short time effect, already evi- 
dent after one iteration as explained above. It illuminates 
another variation on the theme of this paper: on scales 
finer than SZCdV, non-Gaussian statistics may prevail. 
Here, we see a very structured and nonrandom fringe pat- 
tern, which however varies on a scale proportional to h, 
doing no harm to the Schnirelman limit. 



B. Numerical method 



The simple method which we use to find the S'-matrix 
makes use of the expansion of Eq. ^, including up to 70 
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or 100 closed channels along with all the open channels 
(50-500 here) as a basis. Linear equations are set up by 
requiring that '^{x,y) vanish at up to 1000 points {x,y) 
along the right hand wall. The basis functions already 
vanish along y = and y — b, which is a mixed bless- 
ing, since this is also true beyond the right hand wall, 
where this vanishing is unphysical. Without the inclu- 
sion of closed channels the method does not converge. 
The rectangular linear problem ( M by A^', where N' is 
the total number of channels, including the evanescent 
modes, and M is the number of points along the wall set 
to zero) is then solved by singular value decomposition. 
The rationale for inclusion of closed channels is that they 
handle details of the boundary conditions at the disk on 
a scale smaller than a wavelength. The closed channels 
are naturally all taken to have total energy equal to the 
scattering energy, using k'^ = k'^ + ky with k^ pure imag- 
inary and k^ < 0. The values of ky used were given by 
the quantized values in the corridor, however this would 
not be necessary if we included additional points along 
the top and bottom walls near the right hand end, and 
explicitly forced the total wavefunction to vanish there. 

We find that with the restricted basis described above 
the convergence is poor if the disk protrudes too far into 
the billiard. By keeping the center of the disk well to the 
right of the wall, we are able to get stable results for ener- 
gies such that kd < IOtt, where d is the distance the disk 
protrudes. This means that the obstruction can be made 
at least several wavelengths wide in both dimensions, a 
requirement that we must satisfy in order to be in the 
high energy regime. Typically the states we study are in 
the range of the 10,000*^ to 100,000*'' eigenstate of a fixed 
boundary condition billiard; it is possible to go beyond 
the one millionth state for small disks. The range of sta- 
bility of the method may perhaps be greatly extended by 
generalizing the basis to more flexible evanescent modes, 
as discussed above. 

The disk covers a fraction / of the right hand wall. We 
take that fraction to be between 0.04 and 0.28. In anal- 
ogy with the map discussed above, a fraction 1 — / of the 
incoming wave is not scatterered on the first bounce, ap- 
proximately independ ent o f the incoming channel. The 
discussion of Section IV B| holds without modification, 
including the quantum factor of two in the effective 
cross section of the disk, corresponding to diffraction into 
nearby channels. 

The localization of the wavefunction which we now 
have come to expect in channel space ultimately arises 
from the fact that only a small fraction of the incoming 
channel is scattered after each iteration of the S'— matrix 
for small /. The typical scattering channel presents fresh 
amplitude to the disk after each bounce, scattering an- 
other fraction / of the remaining amplitude. The re- 
sulting slow decay out of the initial channel is already 
enough to cause gross anomalies in the wavefunction 
statistics, as compared to RMT. Specifically, this arises 
from the short-time induced Lorentzian envelope in the 
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FIG. 9. The probability distribution for the IPR's is plot- 
ted for various values of /, showing the predicted scaling be- 
havior. 



quasienergy spectrum, as discussed in Sec. IV B 2 



Classically there are also now a finite number of angles 
9, with (a/6) tan — n/m for integer m and n, which 
never hit the disk. For channels corresponding to propa- 
gation near these angles there is a reduction in scattering 
out of the initial channel. These channels are not true 
bouncing ball modes, but near enough to have a strong ef- 
fect on lifetimes. ( "Time" is now the number of iterations 
of the L/- matrix.) Again in complete correspondence with 
the discussion above, the tails of various distributions are 
governed by these near-bouncing ball orbits. 



C. Numerical findings— Sinai billiard 

We consider first the return probability (inverse partic- 
ipation ratio) measures. The scaling relation 7'/(IPR„ = 
x) — fPifx) was predicted in Eq. 40; a plot of fV{fx) 
vs fx for various values of the disk size / is shown in 
Fig. pi confirming this scaling over the whole domain of 
IPR values. We see also from the plot that the typical 
IPR in the Sinai system is « 2//, which for the values 
of / considered is much larger than the RMT-predicted 
value of 3. We also see the expected broad distribution 
of IPR's, with A^— independent width, in contrast to the 
RMT prediction that the spread in the IPR distribution 
should go to zero as I/Vn. 

The tail of the IPR distribution is predicted (Eq. H3) 
to have the power-law behavior 7'(IPR„ = x) ^ 1/fx^. 
This inverse-square behavior was indeed observed, and 
is similar to the same falloff already seen in Fig. g for 
the kicked Sinai maps. The power law tail together with 
the cutoff in the maximum IPR lead to the prediction 
of Eq. ^ namely (IPR„)„ = (P„„) - log TV//. A plot 
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FIG. 11. The distribution of P„^ is plotted for small val- 
ues of P„^ and compared with theory, Eq. ^. In this case 
the bump size is / = 0.072 and the number of channels is 
N = 395. 



of the the dependence of the average IPR on N and / is 
given in Fig. 11^ , where the agreement with Eq. E^ is seen 
to be excellent. As predicted, the mean IPR diverges 
logarithmically away from its ergodic value of 3 in the 
classical limit. 

The distribution of small intensities P„^ should be 
given for our S-matrix by Eq. Efl. This behavior at the 
low end of the P„^ distribution is in very good agreement 
with the theory (Fig. [Tl] ). 

Finally we consider the tails of the intensity and trans- 
port measures. From Eq. pll we expect a cubic fall off in 
the tail of the P„^ intensity distribution: V{Pn( = x) '^ 
1/fx^ for X » 1//. In Fig. ^we display the predicted 
and numerical results, showing good agreement between 
the two. This behavior is controlled by the near-bouncing 
ball dynamics. (We do not discuss again the behavior of 
the intensity distribution intermediate between the head 
and tail; in Fig. H we saw already that the entire distri- 
bution is well-predicted classically.) 



Tail of Pn^ 




FIG. 12. The tail region of the P„^ distribution shows good 
agreement with the predicted cubic power law, for a 27r x 27r 
bilhard, / = 0.23, 226 channels. 




FIG. 13. The tail region of the plot of the P„„' distribution 
shows good agreement with the predicted quartic power law, 
for a 27r X 27r billiard, / = 0.23, 226 channels. 
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The tail of the transport distribution measure V{Pnn') 
is given by Eq. ^ V{Pnn' = x) ^ l/px\ Fig. ^ again 
demonstrates very good agreement with this estimate. 
Notice that the RMT prediction is P„„' — 1 for aU chan- 
nels n ^ n' . 



VII. SIMPLE UNITARY MATRIX MODEL 

The previous examples corresponded to physical sys- 
tems, or maps, which have a direct basis in dynamics. 
Above, we have compared the results for such dynamical 
systems to random matrix theory. However there is a 
variant of the usual random matrix theory, i.e. a modi- 
fied random matrix ensemble, which retains some of the 
gross characteristics of our dynamical systems, while re- 
maining free of any real dynamics. The main idea is to 
retain the tendency to scatter back into the same chan- 
nel (diagonal dominance) while making that portion of 
the amplitude which does scatter do so randomly. This 
gives rise to Lorcntzian envelopes as in the dynamical sys- 
tems, but not with the near-bouncing ball effects, which 
strongly skewed the tails of the intensity, IPR, and trans- 
port distributions discussed above. The S-matrix for the 
collision off the Sinai obstruction in the corridor is the 
key element in the theory of the eigenstates presented 
above. The S-matrix for this process is unitary and sym- 
metric. The fully random matrix ensemble corresponding 
to a symmetric S-matrix is Dyson's circular orthogonal 
ensemble, the COE [g2|. However we wish to modify this 
to a diagonally dominated symmetric unitary matrix U 
which includes the effect of random off-diagonal coupling 
of variable strength. The random component schemat- 
ically represents the scattering off a small object, with 
none of the subsequent dynamical correlations built in. 
We take the form 



U = exp[i{D + R)] 



(72) 



where D is a diagonal matrix with randomly chosen 
quasienergies on the interval [0, 27r), i? is a GOE random 
matrix. Ensembles similar to this have been used previ- 
ously to model spectral statistics intermediate between 
Wigner-Dyson and Poisson |Q. 

For large 7 we approach the COE limit; one iteration 
of the U matrix on a starting vector will then decorrelate 
it completely. That is, the self-overlap becomes 



\{n\U\n)\' 
|(n|n)P 



1 

N 



(73) 



This one step decay corresponds to a pseudoenergy un- 
certainty of 2n, which is just that required to give a uni- 
form spectral density on [0, 27r). For smaller 7 we have 
slower decay, going as 



which leads to a Lorentzian lineshape as in Eq. ^. 

The measures of distributions, tails, etc. above can 
be defined for the Lorentzian envelope (modified COE) 
model as well. The situation here is less rich, since almost 
all states decay with approximately the same rate, unlike 
the special channels (angles) in the Sinai models which 
dictate anomalously slow decay. The IPR is anomalous, 
up by a factor of ~ 7"^ from the RMT prediction due to 
the Lorentzian LDOS envelopes. Transport is similarly 
anomalous, with each channel coupled to only a fraction 
~ 7 of all the other channels at long times {Q ^ I/7). 

Due to the near-bouncing-ball orbits, the average IPR 
in the Sinai systems revealed a localization increasing as 
log TV// as the classical limit N -^ 00 was taken. For the 
unitary matrix model, we have only a single decay rate 
7 (not the distribution caused by the near bouncing ball 
orbits); the resulting average IPR is therefore predicted 
to be independent of N. The decay rate 7 is given by the 
variance of the matrix R through the Golden Rule, 



7 == 2t:{R^)p 



(75) 



with p — N /2-K. In the limit of small 7, the IPR should go 
as 3-2/7, where the factor of three is the COE fluctuation 
factor. The table below shows the mean IPR's averaged 
over all the basis states for two values each of N and 
7 using the modified COE. Excellent agreement is seen 
between the predicted and found IPR's. 



Table of IPR for modified COE 


7|iV 


350 


450 


IPR = 6/7 


0.084 


62.6 


65.9 


IPR = 71.4 


0.188 


32.9 


34.7 


IPR = 31.9 



|(n|C/"|n)p 



exp[— 27m], 



(74) 



The long power-law tails in the intensity, IPR, and 
transport efficiency distributions, present in the real dy- 
namical Sinai system, are similarly absent in the modified 
COE model. 



VIII. CONCLUSION 

Random matrix ensembles possess eigenstates which 
are maximally random, consistent with the symmetry 
constraints which govern the particular ensemble. The 
properties of such eigenstates form the basis of much 
work in quantum chaos theory, and more importantly 
the basis of much theory of nuclei, molecules, and espe- 
cially mesoscopic devices. Is random matrix theory the 
limit to which real classically chaotic systems adhere as 
?i ^ 0? Definitely not. 

The SZCdV theory predicts only coarse grained ergod- 
icity of individual eigenstates in the 7i — > limit, which 
is much weaker than the requirements of random matrix 
ensembles. This gap, between random matrix ensembles 
on the one hand and SZCdV on the other, leaves open 
many questions about the true nature of eigenstates of 
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classically chaotic systems in the ?i ^ limit. We have 
been engaged for some time in the exploration of these 
questions, which address the fluctuations of eigenstates 
on scales that shrink as some positive fractional power of 
Planck's constant (or, more physically, as some negative 
fractional power of the energy). Since such scales become 
infinitesimal as ?i -^ 0, SZCdV has little to say about 
them. Yet they may contain infinitely many wavelengths 
in this limit. The earliest work in this area is scar theory 
P,P| JlJ^Jlq| , which showed that the effects of the least un- 
stable periodic orbits survive the h —^ limit. However 
this is only one possible type of a non-RMT "anomaly" 
in classically chaotic systems. 

Our first investigation beyond scar theory, using the so 
called tilted-wall billiard, examined a very slowly ergodic 
classical system with the expectation that its eigenstates 
would be maximally likely to show non-RMT behavior. 
Indeed the eigenstates did show increasing localization 
on small scales as ?i ^ 0, while still of course obeying the 
SZCdV ergodic theorem. 

In the present study we have switched to the tradi- 
tional paradigm of classical chaos, namely the Sinai bil- 
liard and some close cousins. We have been able to show 
that the eigenstates are ever more strongly localized in a 
certain basis as ?i -^ 0. The basis used is not extraordi- 
nary: essentially it is the usual plane waves (momentum 
space) of scattering theory. We showed that the mean in- 
verse participation ratio in the Sinai-like systems diverges 
logarithmically with increasing energy (or decreasing h) , 
implying that wavefunctions are becoming less ergodic 
at the single-channel scale as the classical limit is ap- 
proached. The situation here is more remarkable than 
in the tilted billiard Q, since in Sinai systems the Lya- 
punov exponent is positive and classical correlations de- 
cay exponentially. A major conclusion of this work is 
that the logarithmically increasing mean IPR in Eq. ^ 
is not due to the bouncing ball states but instead to the 
"near-bouncing ball" channels, whose decay time is large 
compared to the typical decay time 1// but still small 
compared to the Heisenberg time N at which individual 
eigenstates are resolved. 

Another key point is that short time quantum dynam- 
ics and correlation functions have an irreversible effect 
on the localization properties of the eigenstates, as in 
the case of scar theory. 

Undoubtedly there are many more non-RMT effects in 
eigenstates yet be uncovered in other systems, including 
some that could affect important physical properties. 
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